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Abstract 

We argue that the ratio Sj,= ^H/(^He x ^) is a good representation of the local correlation between baryon number 
and strangeness, and therefore is a valuable tool to probe the nature of the dense matter created in high energy heavy- 
ion collision: quark gluon plasma or hadron gas. A multiphase transport model (AMPT) plus a dynamical coalescence 
model is used to elucidate our arguments. We find that AMPT with string melting predicts an increase of 5 3 with 
increasing beam energy, and is consistent with experimental data, while AMPT with only hadronic scattering results 
in a low S 3 throughout the energy range from AGS to RHIC, and fails to describe the experimental data. 
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Data from the Relativistic Heavy Ion Collider (RHIC) 
at Brookhaven National Lab show evidence for partonic 
collectivity and other likely signatures of quark gluon 
plasma (QGP) formation during the early stages of the 
collisions d S S 01 ■ Nevertheless, several important 
questions remain unresolved, such as the beam energy 
where the QGP signatures first appear, and other details 
of the transition between hadronic and deconfined mat- 
ter Quantum Chromodynamics (QCD) predicts a criti- 
cal point separating a first-order phase transition and a 
smooth crossover in the phase diagram of the hot and 
dense QCD matter oila]. It is believed that large fluc- 
tuations in phase space population or large correlation 
length will be one of the experimental signatures of the 
QCD critical point. Investigations of all of these ques- 
tions began at the SPS and the upcoming Beam Energy 
Scan at RHIC |7] will provide an opportunity to study 
them in more detail. Regardless of how difficult or easy 
it will be to uncover specific experimental evidence of a 
critical point, it is a high priority to identify and under- 
stand all the observables that off'er a prospect of discrim- 
inating between hadronic and deconfined matter with 
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good sensitivity. 

In calculations from lattice QCD at high tempera- 
ture, and in models with an ideal quark gas or hadron 
resonance gas, the cross correlations among the con- 
served charges show sensitivity to the confined hadron 
phase or deconfined quark-gluon phase 
Specifically, in lattice QCD |11], the mlw x^ilxl^ the 
strangeness-baryon correlation (x^l) normalized by the 
baryon-baryon correlation Oiff)' approaches unity at 
high temperature in a deconfined phase, and reaches 0.4 
at low temperature in a hadronic phase. The baryon- 
strangeness correlation coefficient Cbs was argued to be 
a robust observable to characterize the nature of the sys- 
tem created in high energy heavy-ion collisions: ideal 
QGP or strongly coupled QGP or hadronic matter 
Although the local baryon-strangeness correlation is a 
sensitive probe of the partonic and hadronic phases as 
predicted by the Lattice QCD ifioll . the proposed exper- 
imental observable (Cbs) is based on global extensive 
quantities 9J. Because it requires a measurement of 
the global baryon number and strangeness in each event, 
an experimental analysis based on Cbs represents a con- 
siderable technical challenge. Further detailed theo- 
retical investigation indicated that a recombination-like 
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hadronization process and hadronic rescattering both 
have the effect of blurring the fluctuation signal Il2lll3ll . 

On the other hand, hypernuclei are clusters of nu- 
cleons and A hyperons lll411 . The production of hy- 
pernuclei happens through a coalescence mechanism 
by the overlapping of the wave functions of protons, 
neutrons and hyperons at the final stage of the colli- 
sions ilSll . This provides a local correlation of baryons 
and strangeness on an event-by-event basis [16J. Specif- 
ically, the deuteron yield is proportional to the baryon 
density while triton (f) and hehum (^He) are a measure 
of baryon correlation llTl llSll . Similarly, hypertriton 
production is related to the primordial A-p phase space 
correlation. The ratio ^3 = ^H/(-^He x ^), which we 
call the Strangeness Population Factor, shows model- 
dependent evidence of sensitivity to the local correlation 
strength between baryon number and strangeness, and 
is demonstrated in this paper to be a promising tool to 
study the onset of deconfinement. The ratio 5 3 is quan- 
titatively a good representation ofx1^/X2 since 5 3 
contains the local strangeness-baryon correlation in the 
numerator and the baryon-baryon correlation in the de- 
nominator UtIi . We expect a prominent enhancement 
of the Strangeness Population Factor in a system that 
passes through a deconfined partonic state, relative to 
what would be observed in a system that always re- 
mained in a hadronic phase |8, 9]. 

In this paper, a multiphase transport model (AMFT) 
is used to study the effects on 5 3 of an existing par- 
tonic phase and the subsequent hadronic scattering. The 
AMPT model is suitable for such studies, since it allows 
us to switch on and off the string melting mechanism 
to simulate a partonic phase when the parton density is 
high at early times, and it also has dynamic transport in 
the early partonic phase and hadronic scattering at the 
late stage. The nuclei and hypernuclei are then produced 
at the final state via Wigner wave-function overlapping 
of their constituent nucleons and hyperon. If the corre- 
lation present at the partonic phase were washed out by 
the hadronic scattering at the later stage, S 3 would have 
been similar, regardless of whether the string melting 
mechanism is on or off in AMPT. In addition, AMPT is 
also used to compute Cbs for a comparison with ^3. 

In a thermally equilibrated system, the yields of 
nuclear clusters via the coalescence mechanism can 
be related to thermodynamic quantities 
However, large fluctuations away from thermal equi- 
librium can result in a locally non-uniform baryon 
and strangeness correlation on an event-by-event ba- 
sis lfl6ll . A dynamical coalescence model has been 
used extensively for describing the production of light 
clusters in heavy-ion collisions [22] at both intermedi- 



ate S El IH and high energies IH, IH H . 
In this model, the clusters are formed in hadron phase- 
space at freeze out. The probability for producing a 
cluster is determined by its Wigner phase-space density 
without taking the binding energies into account. The 
multiplicity of a M-hadron cluster in a heavy-ion colli- 
sion is given by 

Nm = G J dri, c/q,, ■ ■ ■ t/r,„_, t/q,.„ , x 

< 2 ^^(r'l'qi, ■■■'■'«-, 'q,„-,)> (1) 

il>i2>->(M 

In Eq.[T] r,,,- ■ ■,r,„_, and q,_,- ■ ■,q,^^ _ ai-e, respectively, 
the M - 1 relative coordinates and momenta in the M- 
hadron rest frame; pf is the Wigner phase-space density 
of the M-hadron cluster, and (■ ■ ■) denotes the event av- 
eraging. G represents the statistical factor for the clus- 
ter; it is 1/3 for f, ^He fs?) and (A and have the 
same spin as the neutron and triton, respectively). 

To determine the Wigner phase-space densities of 
■'He and ^H, we take their hadron wave functions to be 



those of a spherical harmonic oscillator 1,23^ ,24 125L |2] 

Sin, 



2b^ 



(2) 



and the Wigner phase-space densities are then given by 



exp(-(k2+k2)fc2^ ,(3) 



In Eq. |2] and |3] normal Jacobian coordinates for a 
three-particle system are introduced as in Ref. 11231 124 , 
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. (p,X) and (kp,ki) are the relative coordinates 
and momenta, respectively. The parameter b is deter- 
mined to be 1 .74 fm for ^He S Im S S in and 5 



fm for ^H II33I1 from their rms radii. 



The coordinate and momentum space distributions of 
hadrons (proton, neutron and A) at freeze out are ob- 
tained from AMPT model calculation Q. The AMPT 
model has a good record of agreement with data from 
RHIC [I34I1. including pion-pair correlations llssll and 
flow 113611 . The model has two modes: the default AMPT 
model (version 1.11) involves purely hadronic inter- 
actions only, while the string melting AMPT (version 
2.11) includes a fully partonic stage at the early time 
of the system evolution. Both modes have been used 
in the current analysis in order to distinguish the par- 
tonic and hadronic effect. The overlap Wigner phase- 
space density of the three-hadron cluster, ^H(p,n,A) and 
^He(p,p,n), is then calculated as discussed above, and a 
Monte-Carlo sampling is employed to determine if the 
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Figure 1: (color online) The Wigner phase-space density p for 
from melting AMPT (left panel) and default AMPT (light panel) 
as a function of (A, p) pair momentum. Densities are shown for i/s^ 
= 5 GeV, 17.3 GeV and 200 GeV. The distributions have been normal- 
ized by the number of events at each collision energy. 



cluster is to form a nucleus or not. A nucleus emerges if 
the current sample value p is less than pY^^s^i^y Figure[T] 

depicts the Wigner phase-space density distribution 
as a function of (A, p) pair momentum at various col- 
lision energies in the AMPT model. The coalescence 
probability is larger in melting AMPT than in the de- 
fault AMPT model. The difference increases with colli- 
sion energy. 




Figure 2: (color online) The 53 ratio as a function of beam energy in 
minimum-bias Au-fAu collisions from default AMPT (open circles) 
and melting AMPT (open squares) plus coalescence model calcula- 
tions. The available data from AGS |,331 are plotted for reference. 
The A/p ratios from the model are also plotted. 

Figure |2] shows the S3 results for minimum-bias 
Au-hAu collisions at various beam energies. The defini- 



tion of Strangeness Population Factor (^3 = ^H/(^He X 
j)) incorporates the A/p ratio in order to remove the 
absolute difference in A and p yields as a function of 
beam energy. It is interesting to note that S 3 increases 
with beam energy in a system with partonic interac- 
tions (melting AMPT) while it is almost unchanged in 
a purely hadronic system (default AMPT). The mea- 
surement from AGS 13311 . in spite of large statistical 
uncertainty, gives the value ~ 1/3. The AGS mea- 
surement of 5 4 - ^H/(^He X ^) offers further indi- 



rect support for the lower value of 5 3 at the AGS 03311 . 
A preliminary ^H/^He result for Au-i-Au collisions at 



200 GeV from the STAR collaboration ||37|], in combi- 
nation with the measured A/p ratio from the same ex- 
periment 138, 39, 4^, allows us to infer that the mea- 
sured 5 3 at RHIC is consistent with unity within er- 
rors. These experimental results are consistent with the 
melting AMPT calculations and are in contrast to the 
default AMPT calculations. The data imply that the 
local correlation strength between baryon number and 
strangeness is sensitive to the effective number of de- 
grees of freedom of the system created at RHIC, and 
this number is significantly larger in a system dominated 
by partonic interactions compared with a pure hadronic 
gas. The calculated S 3 from melting and default AMPT 
modes are close at AGS energies and are indistinguish- 
able from the current E864/AGS data. Moving to the 
top SPS energy and beyond, the calculation using melt- 
ing AMPT is more than a factor of two larger than the 
results from default AMPT. It should be noted that if the 
onset of deconfinement takes place at a specific beam 
energy, this may result in a sharper increase of 5 3 than 
the AMPT prediction with string melting scenario. Fur- 
ther experimental efforts are eagerly anticipated, includ- 
ing measurements as part of the RHIC energy scan. 

Furthermore, we investigate the connection between 
our proposed observable 5 3 and the original baryon- 
strangeness correlation coefficient Cbs fl: 



Cbs - -3 



{BS}-{B}{S} 



(4) 



where B and S are the global baryon number and 
strangeness in a given rapidity window in a given event. 
As pointed out in Ref. |fl2ll, a suitable rapidity win- 
dow is important to retain the fluctuation signal. We 
choose the rapidity window of -0.5 < y < 0.5 for the 
present analysis. Figure |3] shows the Cbs in minimum- 
bias Au-i-Au collisions as a function of center-of-mass 
energy from the AMPT model. From top SPS to RHIC 
energy, the Cbs lies between 0.2 and 0.4, and is lower 
than the expected value of unity for an ideal QGP or 
I for a hadron gas [8]. In addition, we find that the 
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Figure 3: (color online) The comparison between S3 and Cbs in 
minimum-bias Au+Au collisions at various beam energies. 



Cbs values from melting AMPT and default AMPT are 
comparable over a wide energy range. As discussed 
in Ref. 1 12], the recombination-like hadronization pro- 
cess itself could be responsible for the disappearance of 
the predicted Cbs deconfinement signal. Detailed study 
indicates that the hadronic rescattering process further 
blurs the signal 1 13(1 . The Cbs increases with an increase 
of the baryon chemical potential //^ 18] at decreasing 
beam energy. The Strangeness Population Factor Sj,, 
on the other hand, increases with beam energy in a sys- 
tem involving partonic interactions, as shown in Fig. |3] 
It carries the potential to reliably resolve the number 
of degrees of freedom of the system created in heavy- 
ion collisions. This suggests that the global baryon- 
strangeness correlation coefficient (Cbs) is less sensi- 
tive to the local baryon-strangeness correlation than the 
Strangeness Population Factor (Sj,) from hypernucleus 
production. Future precise measurements in compari- 
son with our calculations will provide further insight 
into these physics questions that are of central impor- 
tance to relativistic heavy-ion physics. 

In summary, we demonstrate that measurements of 
Strangeness Population Factor S3 are especially sen- 
sitive to the local correlation strength between baryon 
number and strangeness, and can serve as a viable ex- 
perimental signal to search for the onset of deconfine- 
ment in the forthcoming RHIC Beam Energy Scan. 
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